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VERTICAL  AND  ADIABATIC  IONIZATION  ENERGIES  AND 

ELECTRON  AFFINITIES  OF  NEW  SIrC  and  SInO  (ns1-3) 

MOLECULES. 

Alexander  I.  Boldyrev  and  Jack  Simons 
Department  of  Chemistry,  The  University  of  Utah, 

Salt  Lake  City,  Utah  84112,  U.S.A. 

Vyacheslav  G.  Zakrzewski^,  Wolfgang  von  Niessen 
Institut  fur  Physikalische  und  Theoretische  Chemie, 

Technische  Universitat  Braunschweig,  Hans-Sommer  Str.  10, 

D-3300  Braunschweig,  Germany 

Abstract 

Vertical  and  adiabatic  ionization  potentials  as  well  as  electron  affinities 
have  been  calculated  for  SiC,  Si2C,  SisC,  SiO,  Si20  and  SisO  using  five 
different  sophisticated  ab  initio  methods  with  large  basis  sets.  The  geometry 
and  harmonic  frequencies  have  been  calculated  at  the  second-order  Meller- 
Plesset  level.  Results  of  the  calculations  using  all  five  methods  are  in  a  good 
agreement  among  themselves  (±0.3  eV).  The  calculated  vertical  first  IPs  of  SiC, 
Si2C,  SisC  and  SiO  molecules  agree  within  0.2  eV  with  experimental 
appearance  potentials  for  these  species. 


^Present  address:  Department  of  Chemistry,  University  of  New  Mexico, 
Albuquerque,  New  Mexico,  87131,  U.S.A. 
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* 

I.  INTRODUCTION 


Recently,  silicon-carbon  SinC  and  silicon-oxygen  SinO  (ns2,3)  molecular 
clusters  have  been  studied  theoretically^'^  Geometrical  structures  and 
vibrational  frequencies  for  the  ground  electronic  states  have  been  obtained  for 
SiaC  and  SisC  that  agree  well  with  the  available  experimental  data.  In  this 
study,  we  present  ab  initio  calculations  at  several  levels  of  theory  for  vertical 
and  adiabatic  ionization  potentials  and  vertical  and  adiabatic  electron  affinities 
for  SinC  and  SinO  (for  n=1-3).  These  data  may  be  valuable  for  future  mass 
spectrometric  and  spectroscopic  studies  because  several  of  these  molecules 
have  not  been  experimentally  observed.  Hence,  theoretically  predicted 
appearance  potentials  may  be  helpful  for  guiding  experimental  detection.  To 
achieve  as  reliable  results  as  possible,  we  have  used  several  methods  (Meller- 
Plesset  perturbation  theory  (MPn),  quadratic  configuration  interaction 
(QCISD(T)),  conventional  configuration  interaction  (CISD)  and  Green  function 
methods  (the  OVQF  and  ADC(3)  approaches  discussed  below))  and  different 
types  of  atomic  basis  sets  in  performing  the  calculations  whose  results  are 
reported  here. 

II.  COMPUTATIONAL  DETAILS 

The  bond  lengths  of  the  neutral,  cationic  and  anionic  SiC,  SiaC,  SisC, 
SiO,  SiaO  and  SisO  species  were  optimized  by  employing  analytical  gradients^ 
using  the  Gaussian  92  program^  with  polarized  split-valence  basis  sets  (6- 
31 1  +G*)5-7  at  MP2(full)  levels  (UMP2(full)  for  open  shell  systems).  Our  optimal 
geometric  parameters  are  presented  in  Fig.  1  and  our  simulated  IR-spectra 
based  on  the  MP2(full)/6-311+G*  calculated  frequencies  are  drawn  in  Rg.  2. 
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The  MP2(full)/6-31 1  -i-G*  equilibrium  geometries  were  used  to  evaluate  electron 
correlation  corrections  in  the  frozen-core  approximation  by  full  fourth  order® 
Meller-Plesset  perturbation  theory  and  by  the  (U)QCISD(T)  method®  using  6- 
31 1  +G(2df)  basis  sets.  The  UHF  wave  functions  for  open  shell  systems  were 
projected  to  pure  spectroscopic  states  for  which  the  corresponding  results  are 
denoted  PUSCF,  PMP2.  PMP3  and  PMP4.^0 


The  vertical  ionization  energies  and  electron  affinities  were  also 
calculated  by  two  approximations  to  the  many-body  one-particle  Green 
function, namely  the  Outer  Valence  Green  Function  (OVGF)  and  the  third 
order  Algebraic  Diagrammatic  Construction  (ADC(3)). 

In  the  OVGF  method,  all  the  diagrams  up  to  and  including  the  third  order 
terms  which  appear  in  the  expansion  of  the  self-energy  are  included,  and  a 
renormalization  procedure  is  used  for  the  higher  order  corrections.^  The 
renormalization  procedure  includes  three  cases  described  in  details  in  Ref.  15 
and  represents  a  geometric  type  approximation  to  the  self-energy.  Usually  the 
approximation  denoted  case  b)  in  Ref.  15  was  chosen.  Whenever  the  OVGF 
approximation  was  applicable,  alt  three  procedures  yielded  very  similar  results. 

The  ADC(3)  method  takes  all  one-hole,  one-particle,  2  hole  +  1  particle, 
and  2  particle  +  1  hole  configurations  and  their  interactions  into  account. 
Contributions  from  so-called  constant  or  energy  independent  diagrams  were 
evaluated  by  the  iterative  procedure  of  Ref.  16.  The  value  of  the  pole  strengths 
(the  residues  of  the  Green  function  at  each  pole)  are  also  calculated  within  the 
Green  function  method  and  measure  the  validity  of  the  one-particle  picture  of 
ionization  (or  electron  capture).  The  electron  correlation  effects  make  processes 
other  than  the  simplest  Koopmans'  theorem  conization  events  available,  in 
which  case  the  pole  strengths  are  no  longer  equal  to  unity  but  remain  close  to 
unity  when  the  one-particle  picture  of  ionization  (or  attachment)  remains  valid. 


-  -  -  iuia/o? 


When  complete  bteak-down  of  the  one-particle  description  occurs,  the  pole 
strengths  are  much  smaller  than  unity,  and  instead  of  one  photo  electron  line, 
there  appear  a  number  of  additional  lines  known  as  shake-up  (or  shake-off) 
satellites. 

The  SCF  calculations  preceding  the  OVGF  and  ADC(3)  calculations 
employed  the  MOLCAS-2  suite  of  programs^^ .  For  these  calculations,  we  used 
the  very  large  ANO  basis  set  of  Widmark  etal.^^.  For  the  Si  atom,  this  consists 
of  (17s12p5d4f)  elementary  functions  and  [6s5p2d1f]  contracted  functions,  and 
for  the  C  and  O  atoms  of  (14s9p4d3f)  elementary  and  [4s3p2d1f]  contracted 
functions.  The  ANO  basis  set  is  especially  suited  for  the  calculation  of  properties 
where  an  extensive  correlation  treatment  is  necessary  and  an  extended  basis 
set  including  diffuse  functions  is  required.  It  involves  minimal  contraction  loss^^. 
(The  so-called  general  contraction  scheme  is  used  in  these  ANO  basis  sets 
where  all  elementary  functions  enter  each  contracted  function).  These  basis 
sets  are  probably  the  best  for  the  precise  calculation  of  ionization  potentials  and 
electron  affinities. 

The  vertical  and  adiabatic  ionization  energies  and  electron  affinities  were 
also  calculated  for  the  lowest  cationic  and  anionic  states  using  single-reference 
single-and-doubie  excitation  Cl  (CiSD).  The  CISD  energies  were  corrected  by 
the  Davidson  method^o  (CISD(4))  to  account  for  the  quadruple  and  higher 
order  excited  configurations.  These  calculations  were  also  performed  using  the 
MOLCAS-2  program,  where  only  the  core  orbitals  were  kept  frozen  in  all  the 
Green  function  and  CISD  calculations.  The  final  estimates  for  the  vertical  and 
adiabatic  ionization  potentials  and  electron  affinities  are  made  by  averaging  the 
data  at  the  PMP4,  QCISD,  CISD,  OVGF  and-ADC(3)  levels,  and  ±  values  are 
assigned  to  the  range  of  energies  predicted  by  the  various  methods. 


5 


I!'  RESULTS  AND  DISCUSSION 

SIC,  SIC-^  and  SIC'.The  neutral  SiC  molecule  is  known  to  have  a 
triplet  ^Ili  (1 0220^1  ground  electronic  state.2t.22  Our  calculated  bond 

length  Re(Si-C)s1.702  A  and  harmonic  frequency  ve=896  cm't  agree 
reasonably  with  experimental  data:  Re(SI-C)=1.722  A2t  and  964.6  cm*t.2i  An 
experimental  appearance  potentials  has  been  determined  for  SiC:  9.2±0.4  eV 
in  Ref.  23  and  9.0  eV  in  Ref.  24.  We  were  not  able  to  find  in  the  literature  any 
data  on  the  electronic  state,  bond  lengths  or  vibrational  frequencies  of  the 
positive  and  negative  ions  of  SiC. 

From  the  ground  electronic  state  we  can  expect  two  low  energy  cationic 
states:  (1o22o2ijc23ai)  and  2n  {1o22o21jc3)  for  SiC+,  when  an  electron  is 

removed  from  the  highest  a-  or  ic-MO.  We  examined  these  two  states  as  well  as 
the  25;+  (1<j22oM7c^)  state.  A  simple  MO  picture  of  the  valence  MOs  suggests 
the  following:  I0-MO  is  bonding  (sc+ssi) ;  2o-MO  is  antibonding  (sc-ssi) : 

MO  is  bonding  (px.y.c+px.y.si)  and  3o-MO  is  bonding  (pz.c-Pz.Si)-  From  this 
simple  picture  we  expect  that  the  detachment  of  an  electron  either  from  the 
bonding  Itc-MO  (to  give  the  final  state)  or  from  the  bonding  3a-MO  (yielding 
the  2ni  final  state)  should  increase  the  Si-C  bond  length  and  decrease  the 
vibrational  frequency.  In  contrast,  detachment  of  an  electron  from  the 
antibonding  2a-MO  accompanied  by  promotion  of  an  electron  from  3a-MO  into 
the  Itc-MO  (to  give  the  final  state)  should  sharply  decrease  the  Si-C  bond 
length  and  increase  the  vibrational  frequency.  This  trend  is  found  in  our 
calculations  (see  below). 

The  energy  order  of  the  ionic  states  isTiot  easy  to  predict  from  the  simple 
MO  picture.  The  "^X*  state  is  found  to  be  the  lowest  cationic  state.  The  2ni  and 
2X'^  states  are  2.5  eV  and  5.3  eV  higher  at  the  PMP4/6-31 1  +G(2df)  level.  The 
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Optimized  bond  lengths  and  frequencies  are:  1.804  A  and  890  cm'^  for 
1 .675  A  and  952  cm'^  for  2ni;  1 .504  A  and  1443  cm-i  for  25;+  .  We  estimate  the 
first  vertical  IPv=8.9±0.2  eV  for  producing  SiC+  to  be  in  the  range  8.8  eV 
(PMP4),  9.0  eV  (QCISD(T))  and  9.0  eV  (CISD(4)),  which  agrees  well  with  the 
experimental  appearance  potentials  9.210.4  eV23  and  9.0  eV.24  our  first 
adiabatic  I  Pa  is  8.710.2  eV. 

The  experimental  electron  affinity  is  not  known  for  SiC,  therefore  our  data 
may  help  in  the  identification  of  SiC*  anions.  The  two  lowest  anion  electronic 
states  2ni  (1o22a2i7c33o2)  and  2^+  (1<y22a2l  71^30^)  may  be  derived  from  the 
ground  electronic  state  of  neutral  SiC  when  an  electron  is  added  to  the  lowest 
energy  available  MOs.  Again,  from  the  MO-picture,  one  expects  that  bond 
lengths  and  vibration  frequencies  should  decrease  and  increase  respectively, 
for  both  anionic  states. 

The  25;+  state  is  found  to  be  the  lowest  of  SiC*,  however,  the  2ni  state  lies 
only  0.4  eV  (QCISD(T)/6-31 1  +G(2df))  above  the  2£+  ground  state.  The 
optimized  bond  lengths  and  vibrational  frequencies  are:  1.759  A  and  1127  cm*^ 
for  2x+  and  1 .710  A  and  951  cm*^  for  2ni.  Both  of  these  states  are  lower  than 
the  neutral  SiC  molecule,  so  this  molecule  has  two  bound  negative  ion  states. 
The  first  adiabatic  electron  affinity  of  SiC  is  2.25  eV,  which  differs  slightly  from 
the  vertical  EAvs2.32  eV.  Although  EA(SiC)  is  very  high,  it  is  lower  than 
EA(C2)=3.26910.006,25a  EA(C2)=3.27310.00825b  and  close  to 
EA(Si2)=2. 17610.002  eV,26a  EA(Sl2)=2.1991  0.012.26b  jhe  SiC*  ion  also  has 
an  excited  bound  201  state  which  lies  below  the  ground  state  of  SiC  by  1 .83 
eV  and  thus  may  be  found  in  the  gas  phase  experiments.  The  next  state  of 
SiC*  with  1 0^2021 7t^2jci  orbital  occupancy  i&not  electronically  bound  and 
corresponds  to  a  resonance  lying  above  the  energy  of  the  neutral  SiC  by  1 .8 
eV. 
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Si2C,  Si2C+  and  Si2C'.  The  neutral  Si2C  molecule  Is  one  of  the 
major  molecular  species  ( 8%)  observed  in  mass  spectrometric  studies  of  the 
vaporization  of  silicon  carbide  In  the  first  infrared  work,  Weltner  and 
Mcleod^^  proposed  the  tentative  assignment  to  Si2C  of  a  vibration  observed  at 
1187  and  1205  cm*^  in  spectra  of  the  products  of  silicon  carbide  evaporation 
trapped  in  Ar  and  Ne  matrices,  respectively.  Later,  Kafafi  et  al  obtained 
data  which  led  them  to  assign  two  absorptions  at  1 1 88.9  and  658.2  cm'i ,  which 
shifted  to  1 153.7  and  643.3  cm'^  respectively,  and  appeared  to  grow  in  with  the 
same  relative  intensity,  to  the  V3(b2)  antisymmetric  Si-C  stretching  and  vi(ai) 
symmetric  Si-Si  stretching  modes,  respectively,  of  Si2C.  Finally,  a  Fourier 
transform  infrared  study  of  the  vibrational  spectrum  of  Si2C  produced  by 
vaporizing  a  mixture  of  silicon  and  carbon-12  or  carbon-13  has  been  made  by 
Presilla-M^rquez  and  Graham.^^  This  confirmed  a  previously  observed 
vibration  at  1 188.4  cm-f  as  the  V3(b2),  antisymmetric  Si-C  stretching  mode,  and 
resulted  in  the  identification  of  a  new  vibration  at  839.5  cm*f  as  the  vi(ai) 
symmetric  Si-Si  stretching  fundamental.  No  bending  mode  has  been  observed. 
Ab  initio  calculations  performed  by  Grev  and  Schaeferfa,  Bolton  et  aP^  and 
Rittby^c  predicted  a  singlet  angular  C2v  (^Ai)  structure  for  Si2C.  However,  the 
predicted  equilibrium  angle  Si-C-Si  is  very  sensitive  to  the  basis  set^  because 
this  molecule  is  very  floppy  with  a  barrier  for  linearity  of  only  1 .9^®  or  2.1  ^  a 
kcal/mol.  Experimental  appearance  potentials  are  known  for  Si2C  9.210.3  eV^^ 
and  9.1  eV,24  pg  electronic,  geometric  or  vibrational  data  are  available  for 
the  cationic  or  anionic  states  of  Si2C. 

We  optimized  the  geometry  for  the  neutral  Si2C  and  obtained  Re(Si- 
C)=1.707  A,  valence  angle  Si-C-Si=114.7o  and  harmonic  frequencies 
vi(ai)=844  cm-f ,  V2(ai)-129  cm-f  and  V3(b2)=1208  cm*f.  Our  bond  length  and 
valence  angle  reasonably  agree  with  other  CCSO(T)/TZ+2P  (1 .708  A  and 
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1 16.90), ^‘>  and  MBPT2/6-31 1G(2d)  (1.703  A  and  1 19.50)ic  ab  initio  data  if  we 
take  into  account  the  flexible  structure  of  this  species.  Our  harmonic  symmetnc 
and  antisymmetric  stretching  frequencies  vt(ai)=844  cm*^  and  V3(b2)=l208  cm* 
1  also  reasonably  agree  with  experimental  data  vi(ai)s839.5  cm'i,29  and 
V3(b2)= 1188.9  cm-V28 

From  the  1ai2ib222ai22b22lbi23ai2la204ai0  ground  electronic  SCF 
configuration,  we  found  in  OVGF  and  A0C(3)  calculations  that  the  lowest 
energy  cationic  states  are:  2Ai  (1ai2ib222ai22b22 1bi23ai^),  2Bi 
(1ai2ib222ai22b223ai2  ibi^)  and  202  3a^^2bz^)  and  the 

lowest  anionic  states  are:  2A2  (1ai2ib222ai2ibi23ai2ia2^)  and  2Ai 
(1ai2ib222ai22b22lbi23ai24aii). 

In  Koopmans'  approximation,  the  energy  order  of  the  cationic  states  is 
given  as  2Ai<2Bi<2B2.  However,  as  shown  in  Table  I,  at  all  correlated  levels 
and  even  at  the  ASCF  level,  the  order  of  the  cationic  states  is  different: 
2B2<2Ai<2Bi.  Results  of  the  calculations  of  the  vertical  IPs  obtained  at 
correlated  levels  agree  within  0.25  eV  for  all  five  methods:  AMP4,  AQCISD(T), 
ACISD(4),  OVGF  and  ADC(3).  On  the  basis  of  these  date,  we  estimate  the 
vertical  ionization  potentials  to  be:  IP(2B2)=9.1±0.2  eV,  IP(2Ai)=9.3±0.2  eV  and 
IP(2Bi)=9.4±0.2  eV.  The  first  vertical  IP  agrees  well  with  the  experimental 
appearance  potential  for  Si2C  9.2±0.3  eV23  and  9.1  oV.24 

All  of  these  2Ai  ,  2Bi  and  2B2  states  are  well  represented  by  one 
configuration  wave  functions  in  which  the  coefficients  of  the  Hartree-Fock 
determinant  in  CISD  calculations  are  larger  than  0.9  and  the  pole  strengths  in 
the  OVGF  and  ADC(3)  calculations  are  ca.  0.9.  However,  the  higher  ionic  states 
with  ionization  energy  larger  than  12  eV  are  not  well  represented  as  the 
ionization  from  a  single  MO  (see  Table  I).  For  example,  ionization  from  the  next 
2ai-MO  of  Si2C  with  orbital  energy  12.33  eV  is  accompanied  by  strong  many- 


9 


body  effects.  Three  peaks  may  occur  in  the  photoeiectron  spectrum  according  to 
the  A0C(3)  calculations  (see  Table  I)  while  the  configuration  of  the  2ar^  final 
state  remains  strong  (1 1.19  eV)  and  has  a  relative  intensity  of  0.67.  The 
remainder  of  the  intensity  is  shared  by  satellite  lines  among  which  there  is  one 
line  at  13.36  eV  with  a  relative  intensity  of  0.15  (see  Table  I). 

Geometry  optimization  of  the  ^Ai.  ^Bi  and  202  states  of  Si2C+  leads  to  a 
linear  SiCSi  structure.  The  ^At  and  states  collapse  into  the  degenerate 
(1  Og2i  ou22sg2  2au2ijcu3)  state,  and  the  202  state  collapses  into  the  2Xu+ 
(1ag2lou22og2  Ijcu42(yyi).  Por  both  of  these  linear  states,  the  optimal 
geometries  are  presented  in  Fig.  1 ,  and  the  corresponding  adiabatic  ionization 
potentials,  IPa(^nu)  and  IPa(^Iu‘^).  are  equal  to  9.0±0.2  eV  and  9.210.2  eV, 
respectively.  Harmonic  frequencies  have  been  calculated  only  for  2Xu'^  state, 
because  for  the  Hartree-Fock  based  MP2  method,  the  bending  vibrations  are 
not  degenerate.  As  a  result,  even  though  vi(ag)=584  cm*i  and  V2(rtu)='92  cm*i 
appear  reasonable,  the  antisymmetric  stretching  frequency  V3(au)=2405  cm*^  is 
too  high  probably  due  to  symmetry  breaking. 

The  lowest  unoccupied  1a2-  and  4ai-MOs  in  Si2C  have  positive 
energies  at  the  optimal  geometry  of  the  neutral  Si2C  molecule;  therefore,  at  the 
Koopmans'  approximation  these  states  do  not  have  positive  vertical  electron 
affinities.  However,  when  correlation  and  electron  relaxation  are  taken  into 
account,  both  of  these  states  become  bound.  The  vertical  EAy  of  Si2C 
calculated  for  the  2A2  and  2Ai  states  are  0.610.2  eV  and  0.110.2  eV, 
respectively.  Geometry  optimization  of  the  2A2  state  leads  to  a  linear  structure 
and  a  2ng  final  electronic  state,  while  the  finar geometric  structure  of  the  2Ai 
state  is  bent.  Moreover,  in  the  Si2C'  (2Ai)  structure  the  silicon-silicon  distance  is 
very  short  (2.46  A)  and  only  0.13  A  larger  than  a  normal  single  Si-Si  bond 
length.  We  therefore  infer  that  Si2C'  (2Ai)  has  significant  Si-Si  bonding.  The 
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calculated  IR-spectrum  of  SiaC*  (^Ai)  is  presented  in  Fig.  2.  The  adiabatic 
electron  affinities  for  producing  the  two  anion  states,  EAaC^ITg)  and  EAa(^Ai), 
are  0.9±0.2  eV  and  0.7±0.2  eV,  respectively. 

SisC,  SisC-^  and  SisC*.  This  molecule  has  been  detect.^  in  mass 
spectroscopic  experiments  and  an  appearance  potential  of  8.2±0.3  eV  has 
been  reported.^^  jhe  experimental  spectrum  of  SisC  has  been  observed^*^  in  a 
Fourier  transform  infrared  study  of  the  products  of  the  vaporization  of 
carbon/Silicon  mixtures  trapped  in  Ar.  An  ab  initio  stmcture  of  the  SisC  has  been 
determined  by  Rittby^^  at  the  MBPT2/D2P  level,  where  a  rhomboid  singlet 
structure  of  Csv  symmetry  with  carbon-silicon  transannular  bonding  between 
the  two  equivalent  Sib  atoms  has  been  found.  The  ab  initio  vibrational  spectrum 
of  this  stmcture  by  Rittby  is  in  excellent  agreement  with  the  experimental  data  of 
Presilla-Mdrquez  and  Graham.30  No  experimental  or  theoretical  geometries  or 
vibrational  spectra  are  known  for  cationic  or  anionic  states  of  SisC  thus  far. 

We  optimized  the  geometry  and  calculated  harmonic  frequencies  for 
singlet  (1ai2ib222ai23ai2lbi22b224ai2  3b22la202bi0)  rhomboid  neutral 
SisC  (this  geometry  is  shown  in  Fig.  1).  Our  geometry  and  vibrational  spectmm 
(Fig.  2)  for  this  stmcture  agree  well  with  the  MBPT2/DZP  optimized  geometry 
and  corresponding  vibrational  spectmm  by  Rittby. 

We  studied  one  cationic  state  2B2  (1ai2ib222ai23ai2lbi22b224ai2 
3b2^),  which  should  be  the  lowest  positive  ion  according  to  the  vertical  IPs 
calculated  at  the  OVGF  and  ADC(3)  levels.  The  3b2-  MO  does  not  contain 
contributions  from  carbon  AOs,  and  has  antibonding  overlap  between  AOs  from 
the  bridged  silicon  atoms,  however  it  has  a  strong  bonding  interaction  between 
the  terminal  and  bridged  silicons.  Therefore,  detachment  of  an  electron  from  the 
3b2-MO  should  lead  to  a  distortion  In  which  an  increase  of  the  Sit-Sib  bonds 
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might  be  expected.  Indeed,  geometry  optimization  for  the  ^62  state  of  SisC-*- 
leads  to  a  structure  where  the  Sib-Sit  distance  increases  and  the  valence  angle 
SitCSib  increases  from  81°  to  96°  (the  geometries  are  shown  in  Fig.  1  and  the 
simulated  IR-spectra  are  presented  in  Fig.  2). 

We  achieved  agreement  for  both  vertical  IPs  (^62  and  ^Ai)  within  0.3  eV 
for  the  four  different  theoretical  approaches  used:  PMP4,  CISD(4),  OVGF  and 
AOC(3).  From  these  four  methods  our  best  estimates  of  the  two  vertical  IPs  are: 
8.210.2  eV  (202)  and  9.110.2  eV  (2Ai).  The  first  IP  agrees  well  with  the 
experimental  appearance  potential  of  8.210.3  eV.23  The  geometry  of  the  202 
state  differs  significantly  from  the  neutral  molecule,  therefore  the  Franck- 
Condon  factor  for  this  ionization  process  should  be  rather  small  although  the 
adiabatic  correction  (0.4  eV)  to  the  IP  is  not  large  (the  corresponding  adiabatic 
IP(202)  is  7.810.2  eV  (from  PMP4  and  CISD(4))). 

0ased  on  the  MO  energy  ordering  and  results  of  our  OVGF  and  ADC(3) 
calculations,  we  found  two  low  energy  anionic  states:  2A2  (1ai2ib222ai2 
3ai2ibi22b224ai2  3b22la2^)  and  20^  (Iai2ib222ai23ai2ibi22b224ai2 
3b222bi^).  The  orbital  energy  is  negative  for  the  1a2-LUMO  and  positive  for  the 
2bi-MO.  Both  1a2-  and  2bi-MOs  are  pure  rc-type  orbitals.  The  1a2-MO  is 
antibonding  with  respect  to  Sib-Sib  interaction,  but  because  these  two  atoms 
are  situated  far  from  each  other,  occupation  of  this  MO  should  not  lead  to  any 
essential  geometrical  changes.  The  2bi-MO  is  antibonding  with  respect  to  C-Sit 
interaction,  so  elongation  of  the  C-Sit  bond  might  be  expected  in  the 
corresponding  anion.  As  shown  in  Fig.  1 ,  these  types  of  distortions  have  been 
found  upon  geometry  optimization.  The  vertical  electron  affinities  of  these  two 
states,  EAv,  are  estimated  to  bo  1 .310.2  eV  and  0.3 10.2  eV  for  the  2A2  and  20i 
states,  respectively.  Geometry  optimization  does  not  significantly  change  the 
geometry  of  the  SiaC*  (2A2)  state  as  expected;  so  the  adiabatic  EAa=1 .410.2  eV 
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is  very  close  to  the  vertical  electron  affinity.  The  optimal  geometry  of  SisC'  (^Bi) 
is  very  different  from  the  optimal  neutral  geometry  (mainly  because  the  C-Sit 
bond  elongates  by  0.26  A),  so  after  geometry  relaxation  within  C2v  symmetry 
the  EA(2Bi)=0.7  eV  is  higher  than  the  vertical  EA.  However  frequency 
calculations  for  SiaC'  (Cav.  ^Bi)  have  shown  that  this  structure  is  not  a  minimum 
but  a  saddle  point  (with  one  imaginary  frequency).  The  vector  belonging  to  this 
imaginary  frequency  leads  to  an  inplane  distortion.  Beginning  with  a  geometry 
slightly  distorted  from  Cav  (^Bi)  symmetry,  we  refined  the  geometry  optimization 
within  Cs  symmetry.  This  optimization  leads  to  a  Cav  (^Aa)  structure.  With  our 
current  computational  tools  we  are  not  able  to  find  excited  states  for  the 
negative  SiaC*  ion,  but  we  expect  the  existence  of  an  excited  state  on  the  basis 
of  our  calculations  in  the  framework  of  Cav  symmetry. 


SiO,  SiO-^  and  SiO-.  The  SiO  molecule  is  experimentally  known  to 
exist  in  a  singlet  ^2+  ground  state  with  a  bond  length  of  1 .5097  A  and  1241 .56 
cm*^  vibrational  frequency .31  For  SiO+,  Huber  and  Herzberg^i  suggest  1.519  A 
as  the  bond  length  of  the  ground  electronic  state  for  which  the  experimental 
ionization  potential  is  1 1 .6±0.2  eV32  ©r  1 1 .43  eV.33  We  optimized  the  geometry 
of  SiO  at  the  MP2(full)/6-31 1+G*  level  and  found  an  equilibrium  bond  length  of 
1 .536  A  and  a  harmonic  frequency  of  1 183  cm-i. 

The  ionization  energy  of  the  3a>MO  is  smaller  than  that  of  the  1  n-MO  at 
the  Koopmans'  approximation  and  in  the  OVGF  and  ADC(3)  methods.  The 
vertical  ionization  energies  IP(22;+)=:ll.7±0.2  eV  and  IP(2ni)=12.1±0.3  eV  can 
be  considered  as  the  recommended  values  derived  from  the  results  of  the 
present  PMP4,  QCISD(T),  CISD(4),  OVGF  and  ADC(3)  calculations  (see  Table 
I).  The  former  IP(22;+)  is  in  reasonable  agreement  with  experimental  data. 
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Geometry  relaxation  does  not  change  the  lowest  IP  significantly 
IPa(^£'*’)=1 1  -610.2  eV;  the  adiabatic  correction  to  the  second  IP  is  also  small 
IPa(^ni)=11. 710.2  eV.  The  optimized  bond  lengths  and  vibrational  frequencies 
for  these  two  cationic  states  are:  1.534  A  and  1051  cm-i  for  SiO+  (21+),  and 
1 .651  A  and  938  cm*i  for  SiO+  (^rii). 

According  to  our  calculations  (Table  I).  SiO  has  neither  a  positive  vertical 
nor  adiabatic  electron  affinity. 

Si20,  SiaO-*-  and  Si20*.Van  Zee  at  al^^  reported  detection  of  SiSiO 
(X^S)  on  the  basis  of  their  ESR  spectra.  However,  DeKock  et  a|2a  have  shown 
that  the  triplet  SiOSi  linear  structure  is  more  stable  than  the  triplet  SiSiO  by  ca. 
10  kcal/mol  using  CISD/TZ+2P  and  CASSCF/TZ+2P  methods.  Recently 
Boldyrev  and  Simons^^  showed,  at  the  QCISD(T)/6-31 1 -fG(2df)  level,  that 
indeed  the  ground  state  of  this  molecule  has  an  angular  singlet  C2v(^Ai) 
structure  (see  Fig.  1)  while  the  triplet  structures  are  local  minima  and  lie  18-27 
kcal/mol  higher.  The  latter  calculations  are  sophisticated  enough  that  new 
experimental  studies  on  the  structure  of  Si20  molecule  should  be  considered. 

The  assignment  of  the  ionization  spectrum  of  Si20  is  unambiguous  due 
to  the  significant  gap  (about  3.6  eV)  between  the  lowest-energy  ^Ai  (involving 
detachment  of  an  electron  from  4ai-HOMO)  and  the  next  202  state.  The  4ai-MO 
is  G-bonding  and  n-bonding  with  respect  to  Si-Si  and  its  contributions  from  the 
oxygen  atomic  orbitals  are  modest.  Therefore,  after  ionization,  one  expects 
elongation  of  the  Si-Si  bond  length  and  a  relatively  small  change  in  the  Si-0 
bond  lengths.  The  vibrational  vi(ai)  and  V2(ai)  frequencies  in  Si20+  (2Ai) 
which  can  be  attributed  to  the  Si-0  stretch  and  Si-Si  stretch  are  expected  to  be 
approximately  the  same  and  lower,  respectively,  than  the  corresponding  values 
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in  the  neutral  SiaO  molecule.  Results  of  our  calculations  agree  with  this 
qualitative  picture  (see  Rg  1  and  Rg  2). 

The  vertical  IPv(^Ai)  of  SiaO  is  predicted  to  be  7.5±0.2  eV  on  the  basis  of 
the  five  methods  used,  which  is  much  lower  than  the  IP=1 1 .7  eV  of  SiO.  This  is 
because  ionization  in  SiaO  involves  detachment  of  an  electron  from  an  orbital 
residing  mainly  on  the  silicon  atoms,  while  in  SiO,  oxygen  AOs  also  participate 
in  the  orbital.  The  geometric  adiabatic  correction  to  the  first  IP  of  SiaO  is 
relatively  small  (0.3  eV),  and  so  our  estimate  for  the  adiabatic  IPa(^Ai)  is 
7.210.2  eV. 

The  lowest  unoccupied  2bi<MO  has  positive  orbital  energy,  so  the  anion 
is  predicted  to  be  unstable  at  the  Koopmans*  approximation.  This  2bi-MO  is  of 
x-MO  character  with  bonding  interaction  both  with  respect  to  Si-0  and  Si-Si 
overlaps.  However,  Si  atoms  contribute  more  to  this  MO,  so  a  larger  change  is 
expected  in  the  Si-Si  distance  than  in  the  Si-O  bond  lengths  when  this  orbital  is 
occupied.  These  qualitative  predictions  agree  with  the  results  of  our 
calculations  (Rg.  1 )  for  the  Si-Si  bond  length,  but  we  find  the  Si-0  bond  length 
also  increases.  The  vertical  EAv(2Bi)  Is  found  to  be  0.910.2  eV  (see  Table  I) 
and  its  geometrical  adiabatic  correction  is  very  small  (0.1  eV)  so  our  estimate  of 
the  adiabatic  EAv(2Bi)  is  1.010.2  eV.  The  our  harmonic  symmetric  vi(ai)®741 
cm*i  and  V2(ai)s445  cm*^  frequencies  for  Si20*  (^Bi)  look  reasonable,  but  the 
antisymmetric  stretching  frequency  V3(b2)=3917  cm*^  is  overestimated  again, 
probably  due  to  symmetry  breaking. 

SisO,  Si30+  and  SisO'.For  the  cation,  anion,  and  neutral  SiaO,  no 
experimental  data  are  available.  However,  ab  initio  calculations  of  different 
geometric  structures  and  electronic  states  for  neutral  SiaO  have  been 
prevlously2h  performed,  and  a  rhombus  singlet  structure  C2v  (^Ai, 
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1ai22ai2ib223ai2ibi22b224ai25ai22bi23b2°1a20)  found  for  the  ground  state 
(see  Fig.  1 ). 

In  this  work,  we  optimized  geometries  of  the  SisO  cationic  2b i 
(1ai22ai2ib223ai2ibi22b224ai25ai22bi^)  and 
2b224ai22bi25aif)  states,  and  the  anionic  202  (1ai22ai2lb22 
3ai2ibi22b224ai25ai22bi23b2^)  and  2A2  (1ai22ai2ib223ai2ibi22b224ai2 
2bi25ai2ia2^)states,  the  ion  states  having  been  identified  on  the  basis  of 
OVGF  and  ADC(3)  calculations.  The  resulting  geometrical  structures  and 
simulated  IR-spectra  are  given  in  Fig.  1  and  Fig.  2,  respectively.  According  to 
our  calculation,  the  vertical  IPs  and  EAs  (based  on  APMP4,  CISD(4),  OVGF  and 
ADC(3))  are:  IPv(2Bi)=8.1±0.2  eV,  IPv(2Ai)=8.7±0.2  eV,  EAv(2B2)=1  .5±0.2  eV 
and  EAv(2A2)s0.6±0.2  eV  (see  Table  I).  For  all  four  ionic  states,  the  geometries 
do  not  differ  significantly  from  the  equilibrium  geometry  of  the  neutral  SisO 
molecule  (see  Fig.  1).  Therefore,  the  corresponding  adiabatic  energy 
differences  are  IPa(2Bi)*7.8±0.2  eV,  IPa{2Ai)=8.6±0.2  eV,  EAa(2B2)=1.7±0.2 
eV,  and  EAa(2A2)=0.7±0.2  eV. 

IV.  OVERVIEW 

In  this  study  we  used  five  sophisticated  ab  initio  methods  to  calculate 
ionization  potentials  and  electron  affinities:  APMP4,  AQCISD(T),  ACISD(4), 
OVGF  and  ADC(3).  In  the  first  three  methods,  the  ionization  energy  is  calculated 
through  separate  calculations  on  the  neutral  and  ionic  states,  while  two  latter 
methods  compute  the  ionization  energy  directly  as  a  sum  of  electron  relaxation 
and  electron  correlation  corrections  to  the  corresponding  orbital  energy.  Both 
approaches  have  advantages  and  disadvantages.  Indirect  methods  such  as 
AMP4  and  AQCISD(T)  are  based  on  the  unrestricted  Hartree-Fock  reference 
wavefunction  and  therefore  the  resulting  wave  function  does  not  necessarily 
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represent  a  pure  spectroscopic  state.  Direct  methods  have  no  problems  with  the 
purity  of  the  spectroscopic  states  but  cunent  tools  may  be  applied  only  for 
closed-shell  species  whose  ionization  energies  are  to  be  computed.  Therefore 
in  cases  such  as  SiC,  where  the  neutral  and  both  ionic  (positive  and  negative) 
states  are  open  shell,  direct  methods  are  not  applicable.  Another  disadvantage 
of  the  current  direct  methods  is  connected  with  the  calculation  of  adiabatic  IPs 
and  EAs.  However,  promising  rosults35-37  have  recently  appeared  In  the 
literature,  where  geometry  optimization  as  well  as  frequency  calculations  are 
computed  in  direct  methods.  It  seems  that  a  combination  of  the  direct  and 
indirect  methods  is  a  powerful  approach  to  examining  the  ionization  processes. 

Our  ADC(3)  calculations  reveal  for  all  the  species  studied  here  a  break 
down  of  the  one-electron  picture  of  ionization  for  ionization  energies  larger  than 
11-12  eV,  where  many-body  effects  lead  to  satellite  lines  in  the  photoelectron 
spectra.  The  satellite  lines  should  be  more  pronounced  in  the  case  of  SiaO  and 
SisO  than  in  the  SiaC  and  SisC  because  in  SiaO  and  SisO  the  LUMO  and  a  few 
other  unoccupied  MOs  are  lower  in  energy  than  are  the  corresponding  MOs  in 
SiaC  and  SisC. 

Our  calculations  of  the  vertical  IPs  and  EAs  of  SiC,  SiaC,  SiaC,  SiO,  SiaO 
and  SiaO  agree  among  the  five  methods  to  within  0.5  eV.  Our  estimated  first 
vertical  IPs  agree  well  with  the  experimental  appearance  potentials:  8.9±0.2  eV 
vs  9.2±0.4  eV23  and  9.0  eV24  for  SiC;  9.1±0.2  eV  vs  9.210.3  eV23, 9.o24  for 
SiaC;  8.110.2  eV  vs  8.210.3  eV23  for  SiaC;  1 1.710.2  eV  vs  1 1.610.2  eV32  and 
1 1 .4233  for  SiO,  respectively.  We  expect  similar  accuracy  for  the  unknown  SiaO 
and  SiaO  first  IPs. 

The  IP=1 1.7  eV  of  SiO  is  approximately  3  eV  higher  than  the  IP=8.9  eV  of 
SiC.  SiO  has  no  positive  electron  affinity  while  SiC  has  a  large  EA=2.3  eV.  In 
comparison,  Sia  has  IP=7.87  eV,^  and  EA=2.2  eV.26 
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There  are  several  interesting  observations  to  be  made  relative  to  our 
calculated  IPs  and  EAs  of  SinC  and  SinO  and  the  available  experimental  data 
for  pure  silicon  clusters  The  corresponding  carbon  clusters  have  linear 
geometries  and  therefore  are  not  appropriate  species  for  such  comparison. 

When  another  silicon  atom  is  added  to  these  clusters,  things  change; 
IPs7.5  eV  of  Si20  is  lower  than  IP=9.1  eV  of  Si2C,  and  both  of  these  molecules 
have  electron  affinities  near  1  eV.  The  corresponding  Sis  cluster  has  IP=7.9 
eV.41  8.14  eV42  and  EA=2.33  eV.-^a 

Adding  yet  another  Si  atom,  one  finds  that  both  SisO  and  SisC  have 
similar  IPs  (ca.  8  eV)  and  comparable  EAs  (ca.  1.5  eV,  see  Table  I),  and  that  the 
corresponding  Si4  cluster  has  IP=7.6  eV,4i  9.11  eV,^2  7  3  ©744  and  EA=:2.15 
eV.  43  The  IP  and  EA  values  for  SI4  are  close  to  the  corresponding  values  for 
SisO  and  SisC  implies  that  ionization  and  electron  capture  in  SisC  and  SisO 
involves  primarily  the  silicon  part  of  the  cluster. 

As  the  number  of  silicon  atoms  increases  in  SinX  species,  where  X  is 
more  electronegative  atom  than  silicon,  we  expect  that  the  ionization  potentials 
and  electron  affinities  will  become  closer  to  the  corresponding  values  of  pure 
silicon  clusters  Sin-  It  appears  follows  from  our  calculations  that  this  is  the  case 
even  for  rather  small  (n>3)  clusters. 

CONCLUSIONS 

(1)  The  results  of  our  IP  calculations  for  SiC,  Si2C,  SisC  and  SiO  agree 
within  0.3  eV  with  experimental  data  and  among  all  five  of  our  methods  >  PMP4, 
QCISD(T),  CISD(4),  OVGF  and  ADC(3)  to  within  0.5  eV. 

(2)  We  have  made  new  predictions  of  the  vertical  and  adiabatic  IPs  and 
EAs  of  SiC,  Si2C,  SI3C,  SiO,  Si20  and  SisO.  For  the  neutral  and  ionic  states  of 
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these  species,  we  calculated  the  infrared  spectra  to  assist  in  identification  of 
tnese  species  in  the  gas  phase  or  matrix  isolation. 

(3)  We  infer  that  electron  removal  or  electron  attachment  in  SiaC  and 
SiaO  involves  the  silicon  part  of  the  clusters  because  the  corresponding 
ionization  energies  are  close  to  the  corresponding  values  for  Si4-  In  contrast,  for 
Si2C,  SI2O,  SIC  and  SiO,  the  IPs  and  EAs  are  very  different  from  the 
corresponding  numbers  for  Sia  and  Si2- 

(4)  Our  ADC(3)  calculations  reveal  break-down  of  the  one-electron 
picture  of  ionization  at  ionization  energies  higher  than  11-12  eV  where  many- 
body  effects  lead  to  satellite  lines  in  the  photoelectron  spectra.  The  satellite 
lines  should  be  more  pronounced  in  the  case  of  Si20  and  SiaO  than  in  Si2C 
and  SiaC. 
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Figure  Captions 

Fig.  1  Molecular  Structures  and  Geometrical  Parameters  (distances  in  A 
and  angles  in  degrees|at  the  MP2(full)/6'311-i>G*  level  for  SiaC, 
SiaC'*’,  SiaC",  SiaO,  SiaO'*’,  SiaO'.'SiaC,  Si3C+,  SiaC",  SiaO,  SiaO'*’ 
and  SiaO". 

Fig.  2  Calculated  Infrared  Absorbtion  Spectra  for  SiaC,  SiaC-^,  SiaC*, 
SiaO,  SiaO+,  SiaO",  SiaC,  SiaC'*’,  SiaC*,  SiaO,  SiaO'*’  and  SiaO*. 
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